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ABSTRACT 

|Vj ■ Based on the drift-kinetic theory, we develop a model for particle acceleration 

C/^ ■ and transport in solar flares. The model describes the evolution of the parti- 

1^ ■ cle distribution function by means of a numerical simulation of the drift-kinetic 

I . Vlasov equation, which allows us to directly compare simulation results with 

2 ■ observations within an actual parameter range of the solar corona. Using this 

^ . model, we investigate the time evolution of the electron distribution in a flaring 

region. The simulation identifies two dominant mechanisms of electron accelera- 
tion. One is the betatron acceleration at the top of closed loops, which enhances 
the electron velocity perpendicular to the magnetic field line. The other is the in- 
. ertia drift acceleration in open magnetic field lines, which produces antisunward 

^ . electrons. The resulting velocity space distribution significantly deviates from an 

^ . isotropic distribution. The former acceleration can be a generation mechanism 

^ . of electrons that radiate loop-top nonthermal emissions, and the latter be of es- 

caping electrons from the Sun that should be observed by in-situ measurements 



in interplanetary space and resulting radio bursts through plasma instabilities. 
Subject headings: acceleration of particles — plasmas — Sun: flares 



1. Introduction 

Many observations with such as hard X-rays (HXRs), 7-rays, and microwaves have 
revealed that a solar flare is one of the strongest particle accelerator in our solar system. 
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Electrons that are accelerated to several tens of k e V to MeV radiate HX Rs at footpoints 
of flare loops (e.g. 



Sakao 1994: Krucker et al 



above the top of a soft X-ray loop (iMasuda et al 



2003 



1994 



Minoshima et al.ll2009[ ). and sometimes 



19951 ). They are also observed 



via microwaves in the gigahertz band. These emissions have been used to understand the 



prop e rties of accelerated electrons (e.g . . iNakaiima et al. 



2000; Kundu et al. 2001; White et al. 2002 



1983 



Minoshima et al.l 



2008 



Kosugi et al.lll988t ISilva et al. 



). Recently, accelerated 



ions have been studied by the 7 -ray observations of the Reuven Ra maty High Energy Solar 
Spectroscopic Imager {RHESSI; iLin et al.ll2003l : iHurford et al.l 120061 ). 



In addition to these observations of accelerated particles at flare sites, escaping particles 



from the Sun into interplanetary space, the so-called solar energ etic partic 
resulting radio bursts are observed in association with flares (e . g 



es (SEPs), and 



. Lin et al. 


19^ 


^2; 


Lin 


1985 


Krucker et al. 


(2007. 


2009) 



have reported detailed examinations of the relationship of electrons between at the flare site 
and in interplanetary space, by using RHESSI and WIND observations. On th e other hand. 



escap i ng electrons are often observed with no corresponding flare acti vity (e.g.. iPotter et al. 



1980l : iGosling et al.l 120031 : iMorioka et al.l l2007l : Eastwood et al.ll2010l ). This indicates that 



small-scale particle accelerations frequently take place in the corona. 



Based on the magnetic reconnection scenario fIShibata et al.lll995l ). many authors have 
proposed models to explain the particle acceleration in solar flares. The p article acceleration 
in and near the reconnection region has been well discussed. For example, iLityinenkol fll996[) 
has st udied the DC electric field acceleration in reconnecting current sheets. iHoshino et al. 
( I2OOII ) have discussed the nonadiabatic acceleration in the vicinity of the reconnection region 
at which the gyro radius is comparable to the curvature radius of the magne tic field line. 
Fermi acceleration by the contraction of magnetic islands has been studied by iDrake et al. 
( I2OO6I ). In addition to the acceleration near the reconnection region, stochastic accel- 
eration mec h anisrn s in the reconnection downstream region have been proposed by e.g. . 
Miller et al.l (119961 ) (fast-mode magnetohydrodynamic waves) and iTsuneta fc Naitol (Il998l ) 
(oblique shocks), although necessary high-frequency waves for the particle scattering is poorly 
known in the solar circumstance. 

The models greatly contribute to the understanding of how particles are accelerated 
to observed energies. For the complete understanding of many observed features, however, 
we should take into account the particle transport and dissipation processes as well as the 
acceleration in a realistic environment of the solar corona. This is because observed quantities 
do not necessarily reflect the particl e distribution just at the acceleration site, but are the 
convolution of all of these processes ( lAschwandenl Il998l ) . 



Somov fc Kosugil (119971 ) have proposed a model called "collapsing trap", in which the 
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particle acceleration, transport, and dissipation naturally follow from the shrinkage of mag- 
netic loops driven by reconnection. These processes are almost adiabatic in the trap, because 
the particle gyro scale (~ 10~^ s and ~ 10° cm for electrons) is much smaller than the flare 
scale (~ 10^ s and ~ 10^ cm) due to very strong magnetic fields in the corona. The adiabatic 
acceleration may take place everywhere in the trap, which diminishes a serious problem that 



Kane 


1974; 


Miller et al. 


1997) 



The mechanism can 

be regarded as the subsequent acceleration process that occurs after the acceleration near 
the reconnection region, and as the co-existing process with the acceleration in the down- 
stream region. It is noteworthy to note that the model is very similar to the "dipolarization" 
associated with substorms in the terrestrial magnetosphere. 

In the collapsing trap, particles can adiabatically gain ene rgy from convection elec- 
tric fields E = —v x B through VB and inertia drift motions. iKarlicky &: Kosugil (120041 ) 
have studied the electron distribution function in the shrinking loo p, with consi d ering only 
the acceleration due to the VB drift (the betatron acceleration). lAschwandenl (120041 ) has 
analytically c onsidered the e l ectron transport and the resulting impulsive HXR emissions 
in the loop. iGiuliani et al.l (120051 ) have presented a rigorous treatment of the adiabatic 
motion of a single particle in the loop, by numerically solving guiding-center equations of 
motion. iKarlicky fc Bartal ( 120061 ) have performed a guiding-center test particle simulation 
in a cusp-shaped loop obtained from a magnetohydrodynamic (MHD) simulation. However, 
a complete modeling of the evolution of the adiabatic particle distribution in a flare region 
has not been studied yet. 



Qiu et al. 


(2002 


) and 


Asai et al. 


( 


2003. 


2004) 



Liu et al.l ( 120081 ) have observed the correlation between the hardness of the HXR spectrum 
and the strength of convection electric fields along flare ribbons. These observations indicate 
that the electric field greatly influences the distribution of electrons during their travel in the 
corona. The adiabatic model with actual coronal parameters can be tested through a direct 
comparison with these observations, because the model simultaneously describes particle 
acceleration, transport, and dissipation processes in a very wide area of the corona where 
particles are strongly magnetized. 

To understand the particle acceleration, transport, and dissipation processes in solar 
flares, we theoretically study the evolution of the particle distribution based on the numer- 
ical simulation. For this purpose, we newly develop a comprehensive model based on the 
drift-kinetic theory, in which the evolution of the particle distribution function is described 
by means of a numerical simulation of the drift-kinetic Vlasov equation. The simulation can 
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be performed with actual coronal parameters, allowing us to directly compare it with obser- 
vations. In § [2] we describe equations of particles and electromagnetic fields in our model. 
Simulation results are shown in § 121 The simulation identifies two dominant mechanisms of 
electron acceleration; the betatron acceleration at the top of closed loops and the inertia drift 
acceleration in open magnetic field lines. We discuss the results in § HI and then conclude 
the paper in § |5l 



2. Models 

The evolution of a coUisionless plasma is fully described with the Vlasov equation in 
six- dimensional phase space (three dimensions in the configuration and velocity spaces). The 
equation describes the full kinetics of the particle transport along and across magnetic field 
lines, acceleration by electric fields, and the gyromotion, thus can resolve the inertia and 
gyro scales of particles. When the typical scale of interest is much larger than the gyro scale, 
it is rather convenient for focusing on the macroscopic distribution of particles to adopt a 
guiding-center equation (drift kinetics) than the full kinetic equation. In this paper we treat 
the drift-kinetic Vlasov equation to study the time evolution of the particle distribution 
function following a macroscopic change of ambient fields. 



2.1. Basic Equations 



We begin with the relativistic guiding-center equations of motion ( lNorthroplll963l ). 
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where (r, v) is the guiding-center position and velocity of an individual particle, Vd is the 
drift velocity perpendicular to the magnetic field line, e and mo are the charge and rest mass 
of a particle. 
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and Vd is assumed to be much slower than the speed of hght c. In the drift-kinetic theory, 
dependent variables are reduced to {r,u\\,M) by the assumption of gyrotropy. In addition, 
M is no longer a variable because of the conservation of the magnetic moment (eq. ([3])). 

Since observable quantities are the energy and pitch-angle distribution of accelerated 
particles in the case of solar flares, we convert the variables {u\\,M) into (7,/i), where fi = 
[v ■ b) /v is the pitch-angle cosine of a particle. Their time derivatives are written as 
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where we assume E\\ = 0. Using equations ([2]) and (jH), equations (jSj) and are rewritten 
as 
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The first term of the equations represents the betatron (V-B drift) force, and the second term 
is the inertia drift (general form of the curvature drift in the presence of electric fields) force, 
respectively. The third term of equation ([8]) is the magnetic mirror force. The polarization 
drift {dvE/dt in eq. ^) is negligible by the assumption ve/c <^ 1. 



- 6 - 



The time derivative of the position (eq. ([T])) is expressed by the parallel and perpendic- 
ular drift motions to the magnetic field line. The drift motion consists of the Ex B drift ve, 
VB drift, inertia drift, and polarization drift (eq. dl])). We estimate these drift velocities 
in a typical coronal condition. The E x B drift is a motion of fluid described by MHD, 
up to an Alfven velocity ~ 2000 km s~^. The V-B and inertia drifts are ~ (r^/L) w, where 
Tg ~ 10° cm is the gyroradius and L ~ 10^ cm is the characteristic scale of magnetic fields. 
These are much slower than the Alfven velocity. The polarization drift is further slower than 
these drifts. Therefore equation ([T]) can be well approximated into 



dr 



7 



Ve- 



(9) 



Equations ([7]) 



describe the motion of a particle with the guiding-center variables 



We define the particle distribution function /(r, 7,/i) = dN/ {drd'jdn), where dN is the 
number of particles with (7,/^), and drd'jdfi is the volume element in the phase space. Its 
time evolution is described by the follo wing continuity equation of the phase space d ensity 
(the drift-kinetic Vlasov equation; e.g.. IWimmellll984l : iGurnett fc Bhattacharjeell2005l ). 
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2.2. Electromagnetic Fields 



Equation ( iTOl) traces the orbit of the particle distribution function in ambient fields. For 
the determination of the field s, we employ the analytic model of two-dimensional magnetic 
fields of the flare proposed by iLin et al.l (119951 ). Their model is a superposition of potential 
and horizontal fields to impose an X-type neutral line at {x, z) = (0, a), 
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(11) 
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where d is the depth of a dipole moment. The x- and z-axes correspond to the tangential 
and normal directions relative to the solar surface, and the y-axis is along the neutral line. 
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respectively. Note that the above formulae are slightly d ifferent from the original ones, 
because we take the origin at the surface whereas iLin et al.l took it at the neutral hne. 



We introduce an x-position of the footpoint of the magnetic separatrix xj, which is 
magnetically connected to the neutral line. The relation Ay{0,a) = Ay{xf,0) yields 



a 
d 




(14) 



We further introduce an aspect ratio of the magnetic field configuration R = a/xf, which 
yields 

d ~ 2R ■ ^^^^ 
The geometry of the magnetic field is express ed as a function of R, which is treated as one of 
our simulation parameters. Observationally, lAschwanden et al.l (jl996al Jbl) have statistically 
estimated the ratio between the height of the acceleration site and the half length of the 
footpoint distance from a time-of-fiight analysis to HXR data, giving 1.7±0.4. The magnetic 
field configuration with i? = 1.7 is shown in Figure [H 

The fiare evolution can be characterized by the temporal change of the magnetic field 
configuration. Flares frequently sho w a separating niotion of ribbons in wavelengths sensitive 
to th e chromosphere such as Ha (e.g., Qiu et a02002 : Asai et aPboOS . 2004 : Minoshima et al. 
20091 ). Based on the magnetic reconnection mo del, it is interpr e ted a s a chromospheric coun- 
terpart of the rising motion of the neutral line. iTsuneta et al.l (119921 ) reported a growth of a 
cusp-shaped loop, increasing its height with time. To include the evolution in our simulation, 
we change the footpoint position of the magnetic separatrix with time. The time derivative 
(apparent velocity) of the footpoint position is given as a Gaussian profile. 
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From equations (iMl) and (ITS!) . a and R also change with time. 



The temporal change of magnetic fields induces electric fields. Since the electrostatic 
potential can be assumed to be zero in the two-dimensional system, the ^/-component of the 
electric field is determined from the Faraday's law. 
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and other components are zero. The model electric field is a linear function of z, taking zero 
at the surface. This means that the magnetic field lines at the surface are stationary. Since 
the electric field is proportional to dxf/dt, the strength grows till t = tp (we call the period 
as the rising phase), and then is reduced (the declining phase). 



2.3. Simulation Setup 

Using the listed equations, we numerically solve the drift-kinetic Vlasov equation (fTOj) . 
Since we assume two-dimensionality in the configuration space, the equation is reduced 
to a four-dimensional compressive advection equ ation. To solve it, we adopt the R-CIP- 



CSL2 scheme with operator splitting technique ( iNakamura et al.ll200ll ). which simultane- 
ously solves the integrated values of / as well as / itself to keep "subgrid" information and 
to satisfy mass conservation. 

Simulation parameters are d = 1.5 x 10^ cm and R{t = 0) = 1.7. Considering that 
the position of fiare ribbons corresponds to the root o f the magn e tic se paratrix, we set the 



parameters in equation f|T6l) based on observations. iQiu et al.l (120021 ) have observed the 
apparent motion of fiare kernels seen in the Ha line. The velocity reaches several ten to 100 
km s~^ within a time sc ale of ^ 10 s (see Figs. 5 and 6 in their paper). Referring to the first 



spike in Figure 5(b) of iQiu et al.l . we set Vp = 60 km s \ r = 1.77 s, and tp = 5.0 s. The 



simulation runs till t = 10 s. Length (x, z, a, Xf) is normalized to d. 

We limit a calculation domain of the configuration space to the area below the dash- 
dotted line in Figure [H in which the E x B drift velocity calculated with equations fll2p . 
(IT^ . and f[T7]) does not exceed a typical Alfven velocity in the corona (~ 2000 km s~^). In 
other word, our model can not properly describe the electric field around the neutral line, 
where the E x B drift velocity is calculated to be an unrealistic value. Figure [2] shows an 
example of the E x B drift velocity field distribution within the simulation domain. 

Initially, 3 keV isotropic electrons are uniformly distributed in the configuration space. 
We will discuss the assumption of the relatively high initial temperature in § HI The initial 
uniform distribution is reasonable, because the hydrostatic scale height is much longer than 
the size of the simulation domain. At x = Xmax and z = z^i^^ we impose the open boundary 
condition where incoming fiuxes are assumed while outgoing fiuxes are perfectly lost. The 
incoming fiux is fixed to the initial distribution. We also impose another open boundary 
condition at z = 0, where a cold dense plasma is fixed to exclude the possibility of hot 
electrons originating below the surface. The boundary condition at x = is symmetric: 
f{x,z,'j,fi) = /(— X, 2;, 7, — /i). We note that the initial uniform and isotropic distribution 
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is not a steady state solution in our model, because of the different boundary conditions at 
^ = ^max (hot tenuous plasma) and z = (cold dense plasma). 

We use 64 x 96 cells in [x, z) space, and 64 x 64 cells in (7, /i) space. A domain of 7 is 1.01- 
1.20, corresponding to 0.511-102.2 keV electrons. The grids in 7 space are logarithmically 
spaced. 

3. Simulation Results and Interpretation 

Several snapshots of our simulation results are in Figure [31 in which the spatial distri- 
butions of the number of 20 keV and 50 keV omnidirectional electrons at t = 5 s and 10 s 
are presented. In both energies the electron number increases significantly around the top of 
closed loops. The altitude of the peak position of the electron number decreases with time, 
due to the shrinkage of the loops. The distributions in the loops diffuse along the magnetic 
field lines. The spatial diffusion is more prominent for lower energy electrons. Because of 
the boundary condition at z = 0, the number of electrons decreases at the bottom area. 

3.1. Electron Distribution in Closed Loops 

As seen in Figure [3], a strong acceleration takes place around the loop top, at which the 
fastest Ex B drift velocity is observed (see Fig. [2]). The velocity space distribution is of help 
to understand the mechanism around there. Figure H] shows the velocity space distribution 
function taken at the loop top {x, z) = (0, 0.7) at t = 10 s. The horizontal and vertical axes 
correspond to the velocity parallel and perpendicular to the magnetic field line. 

The distribution significantly deviates from the initial, isotropic Maxwell distribution. 
It shows the loss-cone distribution with a loss-cone angle (white lines) defined as a = 
sin~^ B(x, z) /B{xfp,0), where Xi^p is the x-position at the footpoint of a specific field line 
determined from Ay{xFP,0) = Ay{x,z). The distribution also shows a strong enhancement 
of the number of electrons perpendicular to the field line. The perpendicular enhancement 
is energy dependent; the distribution is more anisotropic for higher energy electrons. This 
means that high energy electrons are produced by the acceleration perpendicular to the field 
line, that is, the betatron (V-B drift) acceleration. Since the electron pitch angle becomes 
close to 90 degrees by the betatron acceleration, accelerated electrons are further confined 
and energized at the loop top with time. 

By taking the second-order momentum of the distribution function, we estimate the 
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average energies (temperatures) parallel and perpendicular to the magnetic field line as 



T^{x,z,t) = -J^ d^J ^dfi^f, (19) 

where m| = /i^ (7^ — 1) and = (1 — /i^) (7^ — 1)- The spatial distributions of the parallel 
and perpendicular temperatures at t = 5 s and 10 s are shown in Figure [51 Around the 
loop top, we can see a shght decrease of the parallel temperature and a strong increase of 
the perpendicular temperature relative to an initial temperature 3 keV. The decrease of the 
parallel temperature is due to the loss of electrons inside the loss cone. The perpendicular 
temperature reaches a maximum ~ 15 keV at [x, z) = (0, 0.7) at t = 10 s. The acceleration 
increases the perpendicular temperature by a factor of 5 relative to the initial. 

On the other hand, the rate of the betatron acceleration for electrons with /i = is 
easily obtained from the conservation of the magnetic moment (eq. ([3])), 

= f . (20) 

where 6 = 7 — 1(-C 1) is the kinetic energy. Since the loop at (x, z) = (0, 0.7) at t = 10 s is 
convected from the initial position (x, z) = (0, 0.89) by the E x B drift, a compression ratio 
at the loop top is 

B{t= 10 s) 0.0935 ^ 

B{t = 0s) 0.0156 ' ^ ' 

which is larger than the estimated increase rate of the perpendicular temperature. This is 
because the betatron acceleration less works for electrons with smaller pitch angles. They 
are accelerated rather parallel to the magnetic field line through the inertia drift. 

To consider the acceleration due to the inertia drift in the closed loop, it is convenient to 
introduce the longitudinal adiabatic invariant, u\\L = constant. Here L is the travel length 
of electrons measured along the field line. 



mirror 



^mirror 

,2 



L = 2 ds = 2 ^^1 + {BjB^Ydx, (22) 

^loop— top ^0 

where Xmirror is the x-position at the mirror point. For electrons with the pitch angle close 
enough to the loss-cone angle, Xmirror = xpp- The inertia drift acceleration in the loop is 
expressed as the increase of the parallel energy due to the decrease of the length of the 
shrinking loop. The rate of the inertia drift acceleration is 



/ r \ 2 
£2 / Li 



^2 



(23) 
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By numerically calculating the length of the loop that is initially at [x, z) = (0, 0.89) (the 
same loop as in eq. fl2T]) ). we obtain 



L(t = s) y / 1.076 



L{t = 10 s)y VO-943 



1.3. (24) 



Electrons that experience the inertia drift acceleration escape from the loop top and their 
distribution diffuses along the field line. This yields the perpendicular temperature that is 
not as much as expected from the field compression (eq. ( 1211) ). Nonetheless, the resulting 
velocity space distribution is enhanced perpendicular to the field line (Fig. H]), because the 
rate of the perpendicular acceleration ~ 6 is larger than the parallel ~ 1.3. This is the 
general feature of acceleration with conserving both the magnetic moment and longitudinal 
invariant. 



3.2. Electron Distribution in Open Field lines 

In the previous section, we discussed the acceleration taking place in the closed loop. 
It is interesting to discuss whether the acceleration occurs also in open field lines. Figure [6] 
shows the electron pitch-angle distribution in the open field line near the magnetic separatrix 
at {x,z) = (0.14,0.87) at t = 5 s. The positive (negative) value of the pitch-angle cosine 
corresponds to the sunward (antisunward) direction. 

The pitch-angle distribution is different between sunward and antisunward directions. 
For sunward (/i > 0) the distribution is almost isotropic, because isotropic electrons are 
continuously injected from the upper boundary. For antisunward (/i < 0), on the other 
hand, the distribution deviates from an isotropic distribution. We see lack of electrons with 
/i ~ — 1, also due to the boundary condition. Since we impose the lower boundary condition 
as almost vacuum for high energy electrons, the antisunward electrons are composed by those 
that are injected from the upper boundary (with /i > 0) and then are reflected (change the 
sign of n) by the magnetic mirror force. Injected electrons with /i ~ 1 can not be reflected, 
and are lost at the lower boundary. Therefore electrons with ~ — 1 are not found along 
open field lines. 

We find a slight enhancement of the number of electrons with fi ~ —0.95 and decrease 
with —0.8 < /i < 0, compared to those with > 0. This is an evidence of the parallel 
acceleration to the magnetic field line due to the inertia drift: Electrons are accelerated an- 
tisunward along the curved line by the centrifugal force. This mechanism has been proposed 



for the gene ration of energetic particles in the pl anetary magnetosphere (e .g. JPelcourt et al. 



19961 . 120051 ) . and in quasi-perpendicular shocks flAmano fc Hoshind 120071 ) 
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According to lAmano fc Hoshind (l2007l ). we discuss the resulting pitch- angle distribution 
of this acceleration. In the rest frame, electrons gain energy from convection electric fields 
through the inertia drift. It is convenient to transform the frame into the so-called de 
Hoffmann- Teller frame moving with Ve relative to the rest frame, in which the electric field 
vanishes and hence electrons do not gain energy from convection electric fields. In this frame, 
however, magnetic mirror points move relative to electrons when a magnetic field line has 
curvature. The velocity of the mirror point measured in the de Hoffmann- Teller frame is 

Vrmxxor = -VEtauO, (25) 

where 6 is an angle of the magnetic field line between the departure and mirror points of 
an electron. When Ve is sunward in the rest frame, the mirror point moves antisunward 
in the de Hoffmann- Teller frame. Sunward electrons experience a head-on reflection with 
the mirror points, which increases the electron parallel velocity by 2f mirror and changes the 
direction. 

The changes of the electron pitch angle and speed by the reflection are written as 

/ / A''^ + 2f mirror \ ,.f ^ 

fi = -( , (if /i > 0), (26) 



v' 

(/iV)^ + (l-/iV, (27) 

where variables with (without) a prime are measured after (before) the reflection. The second 
term of equation (|27l) is the square of the perpendicular velocity, which remains constant 
throughout the reflection. The resulting antisunward electron distribution /' is 

/'(l + e',/i') = /.unward(l + e,/i), (if/i'<0), (28) 

where e = 7 — 1 ^ {v/c)^/2. Note that the sunward electron distribution /sunward is indepen- 
dent of yU, /sunward(l + /^) = /sunward(l + c), because it is an almost isotropic Maxwellian. 
With equations and ( E7|) . equation results in 

/'(I + e', /i') = /sunward(l + c' + Afl' y/ e'e^irror + 4emirror), (29) 
where emi„or = (t'mirror/c)V2. 

To determine emirror, we suppose that an electron with fio{> 0) is initially at S, moves 
sunward, and then is reflected at R. The position of the mirror point is determined from 
the conservation of the magnetic moment. 



(30) 
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which depends on the electron initial condition. The velocity of the mirror point (eq. (125|) ) 
depends on the electron initial condition as well as the field configuration, which is evaluated 
as 

l^mirror = -Ve{S) taU [cOS"^ {b{S) ■ b{R)}] . (31) 

We need to determine the mirror point of electrons with different pitch angles, to calculate 
equation 1^ with (|3B. 

Using equations ([7]) - (jH]), we trace the position of electrons backward in time. Electrons 
are set at {x,z) = (0.14,0.87) at t = 5 s, same as employed in Figure [61 Determining the 
mirror points from this test particle simulation, we obtain tan 6' as a function of the electron 
initial pitch angle in Figure [3 The velocity of the mirror point certainly depends on the 
electron initial pitch angle. Using this result, we can calculate the analytic solution of the 
electron distribution (eq. fl2I?l) ). which is presented in Figure |H1 The solution shows an 
enhancement of the number of electrons with fi ~ —0.95, similar to the simulation (Fig. |6]). 
This confirms that the simulation result is certainly interpreted as the above process. 



3.3. Electron Escape from the Sun 

The accelerated electrons due to the inertia drift propagate antisunward along open 
field lines, reach the upper boundary, and finally escape from the simulation domain. It is 
important to consider the number of these escaping electrons and their time evolution. Figure 
[9] shows the time evolution of the spatial distribution of the number of 20 keV antisunward 
electrons. This figure emphasizes the gradient in open field lines (right-half area), to focus 
on the escaping electrons. 

During the rising phase (t = 1 - 5 s) the number of antisunward electrons decreases 
(identified as upward- shifting contours), meaning that the number of escaping electrons 
increases due to the growth of the electric field. The decrease of the electron number is more 
significant at the higher area, because the electric field is a monotonically increasing function 
of z. During the declining phase (t = 5 - 9 s) the number of electrons recovers (identified 
as downward-shifting contours), because the electric field is reduced while constant electron 
fiuxes are continuously injected from the upper boundary. The electron distribution at t = 
9 s is almost same as at t = 1 s, because the electric field distribution is close. 

The inertia drift acceleration takes place in a whole area of open field lines, especially 
near the magnetic separatrix (e.g.. Fig. [6]). We estimate the number of escaping electrons as 

f cscs^pcil, t) = - JJ dxdz j f{'t)df^ - ^^^^ J /(10s)rf/i , (32) 
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where A{t) is the time- varying area of open field fines in tfie simulation domain, and tfie 
integration over tfie configuration space is implemented only in tfie area of open field lines. 
Here we use /(lO s) as a reference, not tfie initial distribution /(O s), because tfie initial 
distribution is not a steady state solution in our model (see § I2.3p . Figure dU] sfiows tfie 
time profile of tfie number of 20 keV escaping electrons. Tfie solid line is obtained from tfie 
usual simulation result. Tfie dasfied line is from tfie result witfi tfie apparent velocity of tfie 
footpoint Vp = 30 km s~^, wfiicfi is a fialf relative to tfie usual. 

Tfie time profile of tfie escaping electrons is similar to tfiat of tfie apparent velocity 
of tfie footpoint, and tfius tfie electric field (eqs. (|T6|) and f|T7|) ). Tfiis is understood as 
follows. A travel time of tfie escaping electrons in open field lines is < 1 s, wfiicfi is mucfi 
faster tfian tfie time scale of tfie temporal cfiange of tfie electric field. Tfie accelerated 
electrons instantaneously reacfi tfie upper boundary and escape witfiout being influenced by 
tfie temporal cfiange of tfie electric field. Tfierefore tfie escaping electrons sfiould directly 
reflect tfie instantaneous configuration of tfie electric field. 

Tfie maximum number of tfie 20 keV escaping electrons is ~ 1% of tfie initial. We 
confirm tfie similar percentage at maximum in 10 - 50 keV. For electrons below 10 keV, 
tfie number is sligfitly small, ~ 0.3%. Comparing tfie solid and dasfied lines, it is found 
tfiat tfie number of tfie escaping electrons is proportional to tfie electric field strengtfi. As 
discussed, tfie rate of tfie inertia drift acceleration is simply proportional to tfie electric field 
(^^mirror oc ve oc E) SO tfiat it is expccted tfiat tfie number of escaping electrons correlates 
witfi tfie electric field. 



Discussion 



Tfie drift-kinetic simulation fias identified two dominant mecfianisms of electron accel- 
eration in solar flares. One is tfie betatron acceleration at tfie top of closed loops, wfiicfi 
enfiances tfie electron velocity perpendicular to tfie magnetic field line. Tfie otfier is tfie 
inertia drift acceleration in open magnetic field lines, wfiicfi produces antisunward electrons. 
Tfie resulting velocity space distribution significantly deviates from an isotropic distribution. 

Tfie betatron acceleration yields more confinement of fiigfier energy electrons at tfie 
loop top, as previously mentioned in tfie collapsing trap model (ISomov &: Kosugil 119971 : 
Karlicky fc Kosugil l2004l : iKarlicky fc Bartal |2006| ) . Tfiey can be a candidat e for radiating 
loop-top n ontfiermal emissions s ucfi as "aboye-tfie-loop-top" HXR sources (IMasuda et al. 



19941. 119951) and microwav e s fe.g..lKundu et al.ll200ll : IWfiite et al.ll2002l : iMelnikov et al 



Huang fc Nakajimall2009l ). iMinosfiima et al.l (120081 ) fiave found tfiat tfie distribution of elec- 



2002 
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trons is enhanced perpendicular to the field line when they are injected into a loop, from 
both observations of the 2003 May 29 fiare and a numerical modeling of the electron trans- 
port with the Fokker-Planck equation. This result can be interpreted that the betatron 
acceleration is the most dominant mechanism for the electron energization in the fiare. 

The betatron acceleration at the loop top is quite reasonable to radiate loop-top non- 
thermal emissions. However, it is widely known that the most of HXR emissions are from 
footpoints of the loop. This indicates a large amount of electrons precipitating there. In the 
model almost all of electrons are trapped in the loop. One of the reasons for this discrep- 
ancy is that the model has not taken into account the pitch-angle scattering. The pitch-angle 
scattering leads electrons into the loss cone, and significantly increases the rate of the pre- 
cipitation. Many authors have argued the effec t of the pitch-angle sca t tering of electrons 
by the Cou l omb collisions with amb ient plasma ( iMelrose fc Brown! Il976l : lAschwanden et al. 



1998, 1999 



Minoshima et al.ll2008l ). Another agency for the scattering is plasma waves. 
In particular, we suggest whistler waves as a possible candidate, because the perpendicu- 
lar temperature anisotropy at the loop top in the simulation (Fig. |5]) is unstable for the 
whistler wave growth. The enhancements of the pitch-angle scattering by the interaction 
with whistler waves might efficiently work on the precipitation and subsequent HXR emis- 
sions at footpoints, similar to the precipitation of radiation belt electrons and diffuse auroral 
emissions at the Earth. We plan further simulations with including the pitch-angle scattering 
as diffusion terms in equation f|TO|) . 

The number of escaping electrons produced by the inertia drift acceleration has been 
estimated up to ~ 1 % of the background (Fig. [10]). The escaping electrons have been ob- 
served in association with fiares, and their number is estimated to ~ 0.1 — 1 % of the observed 



HXR -emitting electrons precipitating into the chromosphere (e.g., lLinlll974l : iKrucker et al. 
20071 ). The observed number of escaping electron is expected to be less than ~ 1 % of the 
background, because the observed number of HXR-emitting electrons should be less than the 
background. Therefore our estimation from the simulation can account for observations. The 
inertia drift acceleration takes place in a whole area of open field lines and contributes this 
high efficiency, as long as curved magnetic field lines keeps moving. It is noted that the ob- 
servation can be understood without a complex phenomenon such as magnetic reconnection 
for the acceleration. 

Because of its ubiquity, we suggest the inertia drift acceleration as a possibl e mechanism 



for producing escaping electrons that are not always associated with fiares (e.g., 


Potter et al. 


1980; 


Goshne et al. 


2003; 


Morioka et al. 


2007; 


Eastwood et al. 


2010 


). 


Morioka et al. 


(2007) 


have reported the so-called "micro type-Ill" radio bursts. Micro type-Ill radio bursts, char- 



acterized by short-lived, continuous, and weak emissions, are thought to be an evidence of the 
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ubiquitous electron acceleration in the solar corona. They have found that the micro type-Ill 
radio bursts are observed when the active regions bordering on coronal holes appear. It may 
suggest that parent electrons are accelerated near the boundary and then escape into inter- 
planetary space. We have shown that antisunward electrons are produced esp ecially near the 



magne tic separatrix, which qualitatively supports the observational results of iMorioka et al. 



(]2007h . 



In § 13.31 we have discussed that the time profile of the escaping electrons by the inertia 
drift acceleration reflects that of the electric field, because the travel time of the electrons is 
faster than the time scale of the temporal change of the electric field. As long as satisfying 
the condition, it is expected that the time scale of the escaping electrons becomes shorter 
following more rapid change of the electric field than we have employed (r = 1.77 s). It may 
explain type-Ill radio bursts that have an elemental time scale shorter than ~ 1 s, while other 
mechanisms may also contribute to the acceleration. The escaping electrons may further be 
accelerated at a higher altitude, for example, the vicinity of the reconnection region and the 
possible turbulent region. The inertia drift acceleration can contribute to supply electrons to 
the higher corona and the interplanetary space where the secondary acceleration may occur. 



Minoshima et al.l (|2009[ ) have reported that the electron distribution in the loop is en- 
hanced parallel to the field line in the 2006 December 13 flare. They have argued that the 
distribution could be formed, if the inertia drift acceleration efficiently takes place. Such 
a distribution has not been formed in our model. In the drift-kinetic model, the betatron 
and inertia drift accelerations take place simultaneously; the former enhances the velocity 
distribution in the perpendicular direction while the latter in the parallel. The resulting 
distribution is determined by combination of these accelerations. We have found in §3. II that 
the rate of the betatron acceleration ~ 6 is more efficient than that of the inertia drift ~ 1.3 
in the closed loop of the model magnetic field. However, it may be expected that the inertia 
drift acceleration overcomes the betatron in different field configurations, and it is possible 
that the resulting electron distribution is enhanced rather parallel to the field line. 

Let us simply consider the relationship between the magnetic field configuration of the 
loop and the distribution of accelerated particles. Suppose that the magnetic field strength at 
the loop top is proportional to its altitude, B oc z~^. When the magnetic field is compressed 
by a factor of Rc, the altitude decreases by a factor of rI^^. We assume that the travel 
length of particles along the field line L is proportional to the altitude. From equations (120|) 
and ( 123|1 . we describe the increase of the perpendicular and parallel energies through the 
compression as 

(33) 
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It is found that particles are accelerated perpendicular rather than parallel to the magnetic 
field line when ,^ > 2. Figure [11] shows the local index ^ of the model magnetic field at 
the apex of loops {x = 0). The index is much larger than 2 at the higher altitude where 
electrons are more accelerated by stronger electric fields (Fig. [3]). This is consistent with the 
perpendicular enhancement of the electron distribution function around the loop top, shown 
in Figure m Meanwhile, we propose from equation (133|1 that it may be possible to enhance the 
particle distribution parallel to the field line when the loop configuration meets < ^ < 2. To 



verify this hypothe sis, we should per 



observations (llnoue et al. 



2008 



' orm s imulations with various field models, e.g., MHD 



simulation re sults flKarlickv fc Bartal 120061) . and re alistic fields extrapolated from Hinode 



Kusano et al.ll2008l ). Toward the understanding of particle 



acceleration in solar fiares through both the observations and simulations, it is critically 
important to study the relationship between a variety of magnetic field configurations and 
the resulting particle phase space distribution. 

Due to the numerical limitation, we have performed the simulation with a relatively high 
initial temperature 3 keV. Although the acceleration efficiency is dependent on the temper- 
ature, it does not mean that the pre-acceleration is necessary for the proposed acceleration 
mechanisms. In fact, we have also performed the simulation with an initial temperature 
1.5 keV, and obtain similar results. This is to be expected, because the adiabatic betatron 
acceleration takes place similarly regardless of the particle energy. For the inertia drift accel- 
eration, the particle velocity should be faster than the mirror point velocity to be refiected. 
Electrons can easily meet this condition. Therefore the acceleration takes place even for 
electrons with the temperature lower than we have employed in this study. 

Let us discuss the effect of the three-dimensionality of ambient fields, for example, the 
shear motion and/or guiding magnetic fields. In the three-dimensional system, the fields may 
evolve so as to generate field-aligned currents. The career of the current should be electrons, 
because of their high mobility. If electrons are inhibited to move along the field line by such as 
the magnetic mirror force and the pressure gradient (acting as a resistivity), parallel electric 
fields will be ge nerated and ac celerate electrons, to carry the current required from the field 
evolution (e.g., lTsuneta]ll995l ). As a result, their distribution will be significantly modified 



from the two-dimensional system, particularly in the parallel direction. The concept of 
the current-voltage relation has been utilize d especially fo r the auroral particle acceleration 
along the magnetic field line at the Earth (jKnightlll973l ). For the electron acceleration in 
the three-dimensional system, it may be essential to consider the parallel electric field and 
the electrostatic potential. 
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5. Conclusion 

Based on the drift-kinetic theory, we have developed a comprehensive model for particle 
acceleration and transport in solar flares. Using this model, we have simulated the time evo- 
lution of the electron distribution in a flaring region. There are two dominant mechanisms 
of electron acceleration. The betatron acceleration takes place at the top of closed loops. 
It can be a generation mechanism of electrons that radiate loop-top nonthermal emissions. 
The phase space distribution of accelerated particles in the loop strongly depends on the 
magnetic field configuration. The inertia drift acceleration also takes place in open magnetic 
field lines, producing escaping electrons from the Sun. The number of escaping electrons esti- 
mated from the simulation can account for the observed number of flare-associated escaping 
electrons. The inertia drift acceleration is caused by a motion of curved field lines, which is 
not necessarily driven by magnetic reconnection. In this sense, we propose the acceleration 
in this study as a mechanism for producing escaping electrons that are not always associated 
with flares. 

We thank to K. Kusano, T. Yokoyama, T. Amano, A. Morioka, and M. J. Aschwanden 
for fruitful discussions and comments. We also thank to anonymous referee for valuable 
comments to improve our manuscript. T. M. is supported by the Grant-in-Aid for Young 
Scientists (B) ^^21740135 and partly by the Grant-in-Aid for Creative Scientific Research 
of MEXT/ Japan, the Basic Study of Space Weather Predication. This work is a part of 
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Fig. 1. — Magnetic field configuration (eq. (fTT]) ) with an aspect ratio R = a/xf = 1.7. The 
white solid lines are the field lines, and the color contour is the field strength. The dashed 
line is the magnetic separatrix. The dash-dotted line is the upper boundary of the simulation 
domain. 




Fig. 2. — E X B drift velocity field distribution calculated with eqs. f|T2|) . f lT3|) . and f fT7|) at 
t = 5 s. The unit of arrows is in km s~^. The black lines are the magnetic field lines, and 
the dashed line is the magnetic separatrix. 
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Fig. 3. — Spatial distribution of the number of 20 keV (left) and 50 keV (right) omnidirec- 
tional electrons at t = 5 s (top) and 10 s (bottom). The white lines are the magnetic field 
lines, and the dashed lines are the magnetic separatrix. 
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Fig. 4. — Velocity space distribution at {x,z) = (0,0.7) at t = 10 s, corresponding to the 
top of a closed loop. The horizontal and vertical axes correspond to the velocity parallel and 
perpendicular to the magnetic field line. The white lines denote the loss-cone angle. 




Fig. 5. — Spatial distribution of the parallel (left) and perpendicular (right) temperatures 
at t = 5 s (top) and 10 s (bottom), calculated from eqs. (fTSll and (fT9|) . The white lines are 
the magnetic field lines, and the dashed lines are the magnetic separatrix. 
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Fig. 6. — Electron pitch-angle distribution in the open field line near the magnetic separatrix 
{x,z) = (0.14,0.87) at t = 5 s. The positive (negative) /i corresponds to the sunward (anti- 
sunward) direction. Different lines represent the electrons with different energies. Number 
density is multiplied by factors for illustration, as shown in the legend on the upper-right 
corner. 
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Fig. 7. — Tangent of an angle of the magnetic field line between the departure and mirror 
points of electrons, as a function of their initial pitch-angle cosine (see eqs. ( 130|1 and ( 13T|) ). 
obtained from a test particle simulation with eqs. ((7]) - (|9]). 
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Fig. 8. — Electron pitch-angle distribution calculated from the analytic solution of eq. ([29 
The format of this Figure is same as of Figure El 
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Fig. 9. — Spatial distribution of the number of 20 keV antisunward electrons at t = 1 s 
(left), 5 s (center), and 9 s (right). The color scale is set so as to emphasize the gradient in 
open field lines (right-half area). The white dashed lines show the contours of the number 
of electrons, with levels of 75%, 80%, 85%, 90% and 95% of the maximum value. 
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Fig. 10. — Time profile of the number of 20 keV escaping electrons (eq. (132!) ). The number is 
normalized to the initial distribution. The solid and dashed lines show the calculation results 
with Vp = 60 and 30 km s~^, respectively. Just after the start of the simulation, electrons 
initially in the loss cone are lost at the footpoint and the number of electrons drastically 
decreases, because the initial distribution is not a steady state solution in the model. This 
appears as the discontinuous jumps at t = 0.5 s. 
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Fig. 11. — Strength of the model magnetic field (dashed line, eqs. (fT2|) - (fT3|) ). and its local 
index ^ (sohd line) at the apex of loops (x = 0) at t = 5 s, as a function of the altitude. It 
is expected from equation ( l33l) that particles are accelerated more perpendicular (parallel) 
to the magnetic field line at which is larger (smaller) than 2 (denoted as a dotted line). 



